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MECHANISTIC  STUDY  OF  THE  EVAPORATION 
OF  LIQUIDS  FROM  POROUS  MICROPARTICLES 


Introduction 


An  understanding  of  the  evaporation  of  hazardous  liquids  from  porous  solids  is  critical 
to  developing  accurate  hazard  assessments  of  liquid  spills.  While  considerable  work  has  been 
conducted  on  the  evaporation  of  liquids  from  solid  beds,  less  attention  has  been  directed  to 
understanding  mass  transfer  from  the  individual  microparticles  constituting  the  solid.  Abramzon 
et  al.  (1973)  conducted  a  theoretical  analysis  on  the  drying  time  of  porous  particles  in  a  fluidized 
bed. 

In  this  study  experimental  evidence  is  presented  that  the  evaporation  of  liquids  from 
microparticles  composed  of  porous  material  can  be  effected  by  the  pore  dimensions  of  the  solid. 
Li  the  present  study  the  evaporation  of  octyl  alcohol  from  porous  silica  aggregate  microparticles 
is  measured.  A  class  of  Syloid  silica,  commercially  produced  by  W.R.  Grace  Chemical 
Company,  is  used  in  this  study  because  of  its  well  defined  porosity  (see  Table  1).  The 
microparticles  in  this  study  vary  in  diameter  from  100  to  150  micrometers  and  the  primary 
particle  size  of  the  silicas  vary  from  3  to  9  micrometers  in  diameter.  Thus  the  experimental 
evaporation  rates  correspond  to  the  evaporation  of  octyl  alcohol  from  an  aggregate  of  silica 
particles. 

The  experimental  system  used  to  study  microparticle  evaporation  is  referred  to  as  the 
single  particle  electrodynamic  levitater  (SPEL).  SPEL  is  an  ultrasensitive  gravimetric  method 
which  uses  levitating  electric  fields  to  detect  the  weight  change  of  individual  particles.  The 
sensitivity  of  the  system  is  such  that  weight  changes  on  the  order  of  tens  of  nanograms  can  be 
detect^  with  this  system.  An  added  feature  of  this  unit  is,  that  by  working  with  microparticles 
100  micrometers  in  diameter,  telemicroscopy  can  be  used  to  monitor  morphology  changes  in  the 
particle,  vis-a  -vis,  wet  and  solid  particles.  This  capability  is  important  when  correlating  the 
liquid  evaporation  rate  to  the  physical  state  of  the  microparticle. 

Phenomenologically,  the  evaporation  of  liquids  from  porous  solids  is  controlled  by 
evaporation  from  the  interstitial  volume  of  the  bulk  solid  (interparticle)  and  evaporation  from 
the  pores  of  the  individual  microparticles  constituting  the  solid  (intraparticle).  In  this  study, 
interstitial  evaporation  is  neglected  and  evaporation  is  assumed  to  be  controlled  by  evaporation 
through  the  pores.  A  shrinking  core  evaporation  model  (Rubel,  1987)  is  developed  and  used  to 
compare  experiment  and  theory.  The  shrinking  core  model  assumes  that  the  liquid  evaporates 
through  an  ever-thickening  porous  shell  which  is  parameterized  by  an  effective  diffusion 
coefficient.  The  effective  diffusion  coefficient  depends  on,  among  other  things,  the  shell  pore 
radius.  Rubel  (1987)  showed  that  for  ammonia  gas  diffusing  through  an  ammonium  phosphate 
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shell,  the  rate  of  gas  diffusion  decreased  as  the  pore  size  decreased.  By  comparing  theory  to 
experiment,  the  effective  diffusion  coefficient  for  the  silica  is  determined.  Using  ^Tieeler's 
(1955)  parallel  pore  model  in  combination  with  capillary  diffusion  theory,  the  effective  pore 
radius  of  the  silica  is  determined  and  compared  to  literature  values. 


Experimentation 


A  schematic  of  the  particle  trap  and  ancillary  equipment  is  shown  in  Figure  1,  The 
particle  trap  consists  of  a  central  electrode  in  the  shape  of  a  hyperboloid  of  revolution  over 
which  an  oscillating  potential  is  applied.  This  electrode  geometry  establishes  an  electric  field 
which  produces  a  net  time-averaged  force  on  a  charged  particle  causing  the  particle  to  execute 
bounded  oscillatory  motion.  Superimposed  on  the  oscillating  electric  field  is  a  static  field 
created  by  two  endcap  electrodes  oriented  such  that  the  static  electric  field  force  counterbalances 
the  gravitational  force  on  the  particle.  A.C.  current  voltages  typically  vary  between  500  and 
1000  volts  and  D.C.  current  voltages  vary  between  0  and  200  volts. 


The  D.C.  voltage  is  adjusted  until  the  trapped  particle  ceases  oscillatory  motion  and 
comes  to  rest  at  the  focal  plane  of  the  telemicroscope.  Under  this  condition,  the  electric  force 
equds  the  gravitational  force  and  the  particle  mass  can  be  determined.  The  telemicroscope  is 
equipped  with  a  scanning  graticule  that  i>ermits  particle  diameter  measurements  to  an  accuracy 
of  ±  5%.  Automatic  positioning  of  the  particle  is  achieved  using  a  servo-mechanism.  The 
servo  uses  a  split  photodiode  detector  to  measure  laser  light  scattered  in  the  forward  direction. 
If  the  particle  migrates  from  the  trap  center,  an  error  signal  is  generated  that  is  converted  to  a 
correction  voltage  that  repositions  the  particle  at  the  center  of  the  trap. 

The  silica-octyl  alcohol  mixture  particles  are  prepared  by  adding  a  sufficient  quantity  of 
alcohol  to  the  powder  silica  until  the  mixture  attains  a  liquid  consistency.  In  some  experiments 
the  silica  was  preheated  for  one  hour  at  an  oven  temperature  of  150°C  and  one  atmosphere 
pressure  to  desorb  preexisting  water  vapor.  However,  evaporation  experiments  showed  no 
difference  between  preheated  and  unheated  silica  samples.  The  mixture  is  stirred  and  allowed 
to  stand  for  a  period  of  one  hour.  After  this  period  of  time,  the  mixture  is  drawn  into  a  pipette 
in  which  a  conducting  wire  is  immersed.  Applying  a  critical  voltage  of  7KV,  a  spray  of  charged 
colloidal  droplets  is  generated.  If  more  than  one  droplet  is  captured  in  the  cell,  a  glass  rod  is 
used  to  remove  all  particles  except  one.  For  the  water  vapor  adsorption  measurements, 
methanol-silica  mixtures  were  prepared  in  the  same  manner  as  the  octyl  alcohol-silica  mixtures! 
Methanol  has  sufficient  volatility  so  that  the  electrospray-generated  particles  would  be  dry  when 
wptured  in  the  cell.  This  flash  evaporation  method  eliminated  the  possibility  of  solvent 
interference  on  water  vapor  adsorption. 

The  gas  environment  of  the  levitated  particle  was  controlled  by  passing  compressed  air 
through  the  levitater  at  a  flow  rate  of  2  cc  per  sec.  This  flow  rate  maintained  the  cell  vapor  free 
of  octyl  alcohol  during  the  evaporation  experiments  and  facilitated  data  analysis.  The  flow  was 
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adjusted  so  that  ventilation  effects  were  minimized.  The  temperature  of  the  cell  was  maintained 
at  25  C  and  all  experiments  were  conducted  at  one  atmosphere  pressure.  For  the  water  vapor 
adsorption  measurements,  the  relative  humidity  of  the  gas  stream  was  controlled  by  mixing  a  dry 
air  flow  with  a  water-saturated  air  flow.  By  varying  the  relative  flow  rate  of  the  two  separate 
flows,  a  precise  control  of  the  gas  dew  point  was  achieved.  The  gas  dew  point  was  monitored 
with  an  EG&G  dew  point  hygrometer. 


Results  and  Discussion 


Water  Vapor  Adsorption 

Initial  experiments  focused  on  the  measurement  of  the  water  isotherm  of  the  silicas. 
Figure  2  shows  the  water  isotherms  for  Syloid  63  and  244  as  measured  by  W.R.  Grace  for  bulk 
powders  and  as  measured  by  the  SPEL  method  for  single  micrometer-sized  particles.  Water 
loading,  the  ratio  of  the  grams  of  water  adsorbed  to  the  grams  of  adsorbent,  is  used  to  quantify 
the  degree  of  water  adsorption.  The  water  loading  is  expressed  in  terms  of  the  levitation  voltage 
according  to  the  relation 


Water  Loading  (gmigm)  =  1  -  /F  , 


(1) 


where  Vq  is  the  levitation  voltage  for  the  dry  adsorbent  particle.  The  water  loading 
measurements  by  SPEL  represent  the  average  of  five  separate  runs,  however,  the  reproducibility 
is  such  that  the  error  bars  are  comparable  to  the  size  of  the  data  points.  Noteworthy  is  the  close 
agreement  between  the  two  distinct  methods  of  measurement,  i.e.,  single  particle  and  bulk 
adsorption.  This  remarkable  finding  suggests  that  SPEL  can  be  used  to  measure  adsorption 
isotherms  of  powders  in  one-tenth  the  time  of  classical  methods,  i.e.,  one  hour  as  compared  to 
one  day.  In  addition,  because  of  the  small  quantities  of  material  required  for  the  measurement 
toxic  material  handling  is  minimized.  ’ 

The  data  in  Figure  2  clearly  shows  that  Syloid  63  adsorbs  more  water  vapor  than  Syloid 
244  for  relative  humidities  less  than  70%.  This  is  primarily  due  to  the  smaUer  pore  size  of 
Syloid  63  which  results  in  pore  condensation  at  the  lower  relative  humidities.  Pore  condensation 
is  the  consequence  of  capillary  filling  which  exhibits  a  concave  meniscus  that  lowers  the  vapor 
pressure  relative  to  a  plane  surface  (Laplace,  1806).  Then  small  pore  condensation  occurs 
before  condensation  onto  the  flat  surfaces  of  the  solid.  .It  is  for  this  reason  that  the  water  loading 
curve  of  Syloid  63  is  elevated  above  that  of  Syloid  244.  ^ 

In  support  of  the  pore  filling  hypothesis,  it  is  noted  that  if  one  monolayer  had  adsorbed 
onto  the  solid  silica,  based  on  the  projected  area  of  a  water  molecule  and  the  specific  surface 
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area  of  the  silicas  as  shown  in  Table  1.,  the  water  loading  for  Syloids  63  and  244  would  be  0.17 
and  0.08,  respectively.  From  Figure  2  it  is  evident  that  more  than  one  monolayer  of  water  has 
adsorbed  onto  Syloid  63.  In  fact  for  the  measured  water  adsorption,  an  equivalent  of  two 
monolayers  of  water  would  have  been  adsorbed  at  60%  relative  humidity.  However,  for  these 
relative  humidities  it  is  unlikely  that  multilayer  formation  is  occurring.  Thus  it  is  postulated  that 
pore  condensation  is  most  probable  in  this  case.  Also  notice  that  for  Syloid  244,  the  monolayer 
calculation  is  in  closer  agreement  with  the  measured  water  loading,  suggesting  that  pore  filling 
is  not  as  active  in  Syloid  244  as  would  be  expected  with  the  larger  pores. 

It  is  interesting  to  note  that  while  the  smaller  pore-sized  Syloid  63  exhibits  greater  water 
vapor  adsorption  at  the  lower  relative  humidities,  the  larger  pore-sized  Syloid  244  shows  a 
greater  saturation  absorption  efficiency  (see  Table  1).  In  fact  the  oil  absorption  efficiency  of 
Syloid  63  is  almost  5  times  greater  than  that  of  Syloid  244.  Thus  as  a  general  rule,  the  efficient 
vapor  adsorbers  are  poor  liquid  absorbers  because  the  greater  internal  surface  area  leads  to  a 
smaller  internal  volume  for  absorption. 

Octyl  Alcohol  Desorption 

Figures  3a-c  show  the  normalized  levitation  voltage  data  for  3  different  evaporating 
mixture  particles  with  the  following  initial  weight  percent  compositions:  a-  38%  Syloid  63,  62% 
octyl  alcohol;  b-  22%  Syloid  244,  78%  octyl  alcohol;  and  c-  25%  Syloid  72,  75%  octyl  alcohol. 
The  initial  composition  is  determined  from  the  initial  and  final  (dry)  levitation  voltages.  Initially 
all  particles  are  trapped  in  the  levitater  as  wet  particles  with  particle  diameters  varying  between 
1(X)  and  160  micrometers.  For  comparison,  the  calculated  levitation  voltage  of  an  evaporating 
droplet  of  octyl  alcohol  with  a  solid  core  and  with  the  same  initial  droplet  diameter  is  shown. 
The  voltage  decay  rate  of  the  pure  droplet  is  calculated  from  the  set  of  equations, 

V  rit  .  '2-nD  Mjd  (P.-  P^ 

;  m  =  - 1 .  ^  (2) 


where  V„  and  m^  are  the  initial  voltage  and  mass,  respectively.  The  mass  decay  rate  is 
controlled  by  continuum  diffusion  and  the  external  octyl  alcohol  partial  pressure  P*  is  assumed 
to  be  zero.  Due  to  the  small  latent  heat  of  evaporation  and  evaporation  rates,  isothermal 
conditions  prevail.  It  is  evident  from  the  calculations  that  the  mixture  particles  evaporate  like 
pure  droplets  initially.  However,  the  mixture  particle  evaporation  rate  diverges  from  the  pure 
droplet  evaporation  rate  during  the  evaporation  process.  From  telemicroscopic  observations,  it 
is  observed  that  the  divergence  in  the  particle  evaporation  rates  is  associated  with  the  particle 
transforming  from  the  wet  to  the  dry  state. 

To  model  the  change  in  the  particle  evaporation  rate,  a  shrinking  core  evaporation  model 
is  introduced  (Rubel,  1987).  The  shrinking  core  model  accounts  for  the  diffusion  of  the  octyl 
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alcohol  vapor  through  the  porous  silica  shell  as  the  octyl  alcohol  liquid  core  shrinks.  The 
shrinking  core  model  assumes  that  the  evaporating  liquid  forms  a  spherical  core  concentric  with 
the  solid  particle  surface.  Evaporation  of  the  liquid  core  is  controlled  by  the  diffusion  of  alcohol 
vapor  through  a  porous  shell  that  has  been  depleted  of  Uquid  by  octyl  alcohol  evaporation.  As 
the  liquid  core  continues  to  evaporate,  the  shell  thickness  increases.  The  diffusion  of  gas 
through  the  shell  depends  on  a  variety  of  parameters  including  shell  thickness  and  shell  pore 
radius.  We  begin  by  developing  the  vapor  concentration  field  within  the  porous  sheU.  If  one 
draws  an  arbitrary  sphere  of  radius  r,  the  surface  of  which  lies  within  the  porous  shell,  the  flux 
through  this  sphere  surface  is  given  by 


47tr2Z)  ^  s  N 
*  dr 


(3) 


where  N  is  a  constant  for  steady  state  conditions.  The  effective  diffusion  coefficient  of  the 
porous  shell  and  the  octyl  alcohol  vapor  concentration  inside  the  shell  Cj  are  shown  explicitly. 
Integrating  Eq.3  from  R-t  to  R,  where  t  is  the  thickness  of  the  porous  shell,  and  using  the 
boundary  conditions  Ci(R-T)=c°,  where  c°  is  the  saturation  vapor  concentration,  N  can  be 
expressed  as 


N  = 


4nD^R  (R-t) 

T 


(4) 


Likewise  by  integrating  Eq.3  from  R-t  to  r,  the  radial  dependence  of  the  concentration 
field  is  expressed  as 


c,(r)  -  c" 


N  f 

4izD^  [  R-x 


(5) 


Substituting  Eq.4  into  Eq.5,  Cj  is  expressed  explicitly  as 
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c,(r)  -  c"  =  (  c,(R)  -  c"  )  . 


To  determine  Ci(R),  flux  continuity  across  the  R-sphere  surface  interface  is  invoked,  i.e., 


4nR2i)  I  R  =  4t:R^D  I  R  , 

dr  ^  dr 


where  the  external  vapor  concentration  fleld  satisfles  the  conditions 


=  cj^oo)  as  r^  oo  . 


Deriving  the  solution  for  c,,  as 


-  cj«>)  )  -  +  cj~) 


and  substituting  into  Eq.7  along  with  Eq.6,  Cj  (R)  is  found  to  be 


c.  (R)  = 


X  -  +  D 


where  it  is  assumed  that  c„(<»)=0.  Finally.  substituUng  Eq.8  into  Eq.6,  the  internal 
concentration  field  is  expressed  as 


D^(R-x)  + 
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The  evaporation  rate  is  determined  by  substituting  Eq.9  into  Eq.3,  i.e., 


—  •> 

V,  =  47rv^2Z).-^ 


47iv^  iR-x)Dp^c^ 
DJiR-x)  +  D^x 


(10) 


Note  that  if  t— 0,  Eq.lO  reverts  to  the  expression  for  pure  droplet  evaporation  as  expected.  To 
integrate  10,  r  must  be  expressed  as  a  function  of  time  or  as  a  function  of  the  liquid  volume 
V|.  Focusing  on  the  shrinking  core  of  radius  R-t,  the  volume  of  the  core  is  expressed  in  terms 
of  the  solid  and  liquid  volumes  as 


V,  +  V,  = 


4ti 


(R-c? 


(11) 


where  v,  and  v,  are  the  volumes  of  liquid  and  solid,  respectively,  within  the  shrinking  core.  The 
liquid  is  assumed  to  fill  the  interstitial  volume  and  the  pores  of  the  primary  particles.  Rewritting 
Eq.ll  ^d  using  the  definition  v/v,  =  l/ep„  where  e  and  p,  are  the  total  free  volume  (cc/gm) 
and  solid  density  (gm/cc),  respectively,  the  shell  thickness  is  written  as 


T 


=  R  - 


(3(1  +  ep;  y/3 


(12) 


Substituting  Eq.l2  into  Eq.lO,  the  liquid  volume  evaporation  rate  is  given  by 


= 


4ir/f  Vi  c° 


\l/3 


3(1+  ep;  Dlj 


.-1/3 


(13) 


Using  separation  of  variables,  Eq.l3  is  readily  integrated  to  give 
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(14) 


Wn 


(yft)  -  v/O))  +  ^  v/O)^  (v/O^  -  v/0)2/3)  = 


where 


v/0) 


47t/?^€p^ 

3(1+  €p^)  ' 


To  compare  the  levitation  voltage  data  to  model  predictions,  the  liquid  volume  v,  must 
be  related  to  the  levitation  voltage  V.  Because  the  relative  levitation  voltage  is  equal  to  the 
relative  particle  mass,  i.e., 


V(t)  ^  +  m, 

V(P)  +  m/O) 


then  using  Eqs.  1 1 
relation 


and  12,  the  relative  levitation  voltage  is  related  to  the  liquid  volume  by  the 


Vjt) 

F(0) 


\  ^  3(1+  ep^P/V/^ 
.  ^P/  ; 


(15) 


pus  levitation  voltages  are  calculated  by  first  solving  Eq.l4  and  then  substituting  the  results 
into  Eq.l5.  The  following  parameter  values  are  used  in  the  calculation;  D  =0.06  cmVsec* 

^>=2.29  gm/cc;  p, =0.8270  gm/cc;  and  v,=157.47  cc/mole.  The  gas 
p  ase  diffusion  coefficient  Dg  is  computed  using  the  molecular  theory  of  gases  for  rigid  sphere 
collisions  (Rubel,  1984).  The  particle  radius  R  at  the  onset  of  shrinking  core  evaporatkm  is 
determmed  by  noting  the  time  at  which  the  particle  transforms  from  the  wet  to  dry  state  and 
usmg  Ae  pure  droplet  flux  rate  (Eq.2)  to  evaluate  the  corresponding  particle  radius.  The 
remammg  parameters  to  be  defined  are  the  free  volume  e  and  the  shell  diffusion  coefficient  D  . 
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The  free  volume  is  determined  from  the  levitation  voltage  of  the  octyl  alcohol-saturated 
particle  V,  (voltage  at  particle  darkening)  and  the  levitation  voltage  of  the  dry  particle  (v,=0) 
using  Eq.l5.  Setting  v,=0,  the  free  volume  is  given  by 


€ 


(16) 


where  Vd  is  the  levitation  voltage  of  the  dry  particle.  From  Figs.  3a-c,  the  free  volumes  are 
found  to  be  3.26,  1.96,  and  1.09  cc/gm  for  Syloids  244,  72,  and  63,  respectively.  These  results 
can  be  compared  to  the  oil  absorption  data  reported  by  W.R.  Grace  (Table  1)  by  converting  the 
free  volume  to  an  oil  weight  equivalent.  Because  the  oil  absorption  data  is  referenced  to  a 
dibutylphthalate  liquid  with  a  bulk  density  of  1.047  gm/cc,  the  free  volume  is  converted  to  an 
oil  weight  equivalent  by  multiplying  the  free  volume  by  the  density  factor  1.047.  The  oil  weight 
equivalent  based  on  the  experimental  data  become  3.41,  2.05,  and  1 . 14  for  Syloids  244,  72,  and 
63,  respectively.  These  oil  absorption  values  compare  favorably  with  the  data  in  Table  1,  i.e,, 
2.80,  2.20,  0.60. 

Once  the  free  volume  is  determined,  the  only  remaining  adjustable  parameter  is  the 
effective  shell  diffusion  coefficient  which  is  evaluated  by  a  "best  fit"  between  data  and 
shrinking  core  model  calculations.  Figures  4a-c  show  model  calculations  for  three  different 
values  of  D*  for  each  of  the  three  mixture  silica  particles.  The  value  DE2  represents  the  "best 
fit"  which  is  determined  by  calculating  the  correlation  coefficient  for  the  data-model  pair.  The 
value  DEI  represents  the  case  where  the  shell  diffusion  coefficient  equals  the  gas  phase 
coefficient,  i.e.,  unhindered  diffusion.  The  calculation  with  DE3  shows  the  sensitivity  of  the 
model  on  the  value  of  the  effective  diffusion  coefficient.  Fig. 4a  shows  the  comparative  fit  for 
the  evaporation  of  the  Syloid  72  mixture  particle.  The  normalized  levitation  voltage  is 
referenced  to  the  saturation  voltage  V,.  The  best  fit  between  data  and  model  is  obtained  with 
DE2=0.019  cmVsec  with  a  correlation  coefficient  of  0.994.  A  similar  analysis  is  applied  to  the 
Syloid  244  evaporation  data  as  shown  in  Fig.4b.  Best  fit  is  found  for  an  effective  diffusion 
coefficient  equal  to  0.025  cmVsec  with  a  correlation  coefficient  of  0.997.  For  the  Syloid  63 
mixture  particle,  the  "best  fit"  value  for  the  diffusion  coefficient  is  0.015  cmVsec. 

Thus  the  ranked  order  for  the  shell  diffusion  coefficient  in  descending  order  is  0.025 
cmVsec,  0.019  cm^/sec,  0.015  cmVsec  for  Syloids  244,  72,  and  63,  respectively.  Qualtitatively 
the  decrease  in  the  diffusion  coefficient  with  decreasing  Syloid  number  is  consistent  with  the 
porosity  data  shown  in  Table  1.  Specifically,  from  Table  1  it  is  evident  that  the  pore  radius, 
which  is  proportional  to  the  ratio  of  the  pore  volume  to  the  surface  area,  decreases  with 
decreasing  Syloid  number.  By  Knudsen  flow  considerations  the  shell  diffusion  coefficient  is 
proportional  to  the  pore  radius  and  thus  according  to  Table  1,  the  pore  diffusion  coefficient 
should  decrease  with  decreasing  Syloid  number  which  is  consistent  with  our  results. 

Quantitatively,  the  pore  radius  can  be  calculated  from  the  experimental  shell  diffusion 
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coefficient  using  Wheeler’s  (1955)  parallel  pore  model  in  combination  with  Knudsen  capillary 
flow  considerations.  According  to  Wheeler’s  parallel  pore  model,  the  shell  diffusion  coefficient 
is  related  to  the  single  pore  diffusion  coefficient  Dp  by  the  relation 


(17) 


where  ij  is  the  surface  porosity  of  the  shell  and  k  is  the  tortuosity  factor  describing  the  deviation 
of  the  pores  from  ideal  cylinders.  The  surface  porosity  of  the  shell  is  defined  as  the  ratio  of  the 
pore  surface  area  to  the  shell  surface  area  and  ranges  in  value  from  0  to  1.  Without  supporting 
data,  the  surface  porosity  is  set  arbitrarily  at  an  average  value  of  0.5  for  all  Syloids.  In  the 
present  calculation  it  is  assumed  that  cylindrical  geometry  exists  and  that  k=1.  If  the  capillary 
radius  is  much  less  than  the  mean  free  path  of  the  gas,  as  is  the  case  for  the  silicas,  then 
diffusion  through  the  capillary  is  described  by  Knudsen  flow.  The  individual  pore  diffusion 
coefficient  for  a  cylindrical  pore  is 


4r 

p  3 


(18) 


where  r^  is  the  pore  radius  and  M„  is  the  molecular  weight  of  the  diffusing  gas.  Then  the  pore 
radius  can  be  determined  from  Eqs.  17  and  18  from  knowledge  of  the  effective  diffusion 
coefficient  and  the  surface  porosity.  Performing  these  calculations  the  pore  radius  is  determined 
to  be  374,  284,  and  224  A  for  Syloids  244,  72,  and  63,  respectively.  Average  pore  sizes  can 
be  calculated  from  Table  1  using  the  pore  radius  expression  given  by  Alishusky  and  Fournier 
(1990)  re=  4Vp/Ap ,  where  Vp  and  Ap  are  the  pore  volume  and  surface  area,  respectively.  Using 
the  values  in  Table  1,  the  literature  values  for  the  pore  radii  are  181,  47,  and  24  A  for  Syloids 
244,  72,  and  63,  respectively. 

It  is  evident  that  the  present  study  overestimates  the  pore  radii  of  the  Syloids  when 
compar^  to  the  porosity  data  of  W.R.  Grace.  One  possible  explanation  is  the  assumption  that 
interstitial  volume  evaporation  can  be  neglected  is  invalid.  Because  the  evaporation  flux  from 
the  interstitial  volume  is  greater  than  that  from  the  pore  volume,  the  presence  of  interstitial 
volume  will  increase  the  best  fit  value  for  the  effective  diffusion  coefficient.  Consequently,  the 
value  of  the  pore  radius  is  increased  due  to  presence  of  interstitial  evaporation.  Further  evidence 
for  the  contribution  of  interstitial  evaporation  is  found  in  the  primary  particle  diameter  data  given 
in  Table  1.  Because  the  aggregate  with  the  largest  primary  particles  possesses  the  largest 
interstitial  volume,  it  is  expected  that  Syloid  63  would  show  the  greatest  disparity  in  the  pore 
radii  and  Syloid  244  would  show  the  smallest.  Indeed  this  is  the  case  as  disparties  in  the  pore 
radii  range  from  a  factor  of  2  for  Syloid  244  to  factor  of  10  for  Syloid  63. 
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Future  efforts  will  be  directed  to  developing  a  composite  model  which  can  account  for 
combined  interstitial  and  pore  volume  evaporation.  This  model  will  be  applied  to  the  present 
data  to  ascertain  a  more  accurate  estimate  of  the  pore  radii  of  the  solids.  Another  approach  is 
to  measure  the  evaporation  flux  from  microparticles  composed  of  single  primary  particles.  This 
would  eliminate  the  presence  of  interstitial  evaporation  and  would  be  amenable  for  analysis  by 
the  shrinking  core  evaporation  model. 


Conclusions 

The  evaporation  of  octyl  alcohol  from  single  microparticles  composed  of  an  aggregate 
of  primary  silica  particles  is  measured  using  single  particle  electrical  levitation.  It  is 
demonstrated  that  the  evaporation  flux  is  retarded  by  the  porous  silica  particles  and  the 
evaporation  retardation  increases  with  decreasing  silica  pore  size.  A  shrinking  core  evaporation 
model  is  developed  that  accounts  for  pore  evaporation  and  effective  diffusion  coefficients  are 
determined  by  comparing  model  results  to  data.  Using  Knudsen  flow  considerations,  pore  radii 
are  evaluated  from  the  diffusion  coefficients  and  compared  to  literature  data.  It  is  found  that 
the  pore  radii  determined  in  the  present  study  exceed  literature  values  suggesting  that  interstitial 
evaporation  (evaporation  from  the  volume  between  the  particles)  is  contributing  to  the 
microparticle  evaporation  rate.  It  is  also  demonstrated  that  single  particle  levitation  can  be  used 
to  determine  the  liquid  absorption  capacity  of  solid  particles  rapidly  and  accurately. 
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Table  1. 


Typical  Physical  and  Chemical  Characteristics  of  Syloid 
(W.R.Grace  Chemical  Co.) 


Silicas 


Syloid 

63 

Si02(%) 

99.4 

Avg. Particle  Size 

(jim) 

9 

Surface  Area 
(mVgm) 

675 

Oil  Absorption 
ab/lOOlb) 

60 

Bulk  Density  (Ib/ft^) 

29 

Soluble  Salts  (%) 

0.5 

Pore  Volume 
(cc/gm) 

0.4 

72 

244 

99.4 

99.4 

4 

3 

340 

310 

220 

280 

11 

8 

0.3 

0.2 

0.4 

1.4 

18 


cc 
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Figure  1.  Schematic  of  Single  Particle  Electrodynamic  Levitation  (SPEL)  Apparatus 
and  Ancillary  Gas  Conditioning  and  Particle  Detection  Devices. 
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Figure  2.  Comparison  of  Syloid  Silica  Water  Isotherms  as  Measured  by  W.R.  Grace 
for  Bulk  Samples  and  as  Measured  by  SPEL  for  Single  Particles. 
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Figure  3b.  Evaporation  Data  for  a  Mixture  Particle  Initially  Composed  of  22  %  Syloid 
244  and  78%  Octyl  Alcohol  and  Calculated  Evaporation  Rate  of  an  Octyl  Alcohol 
Droplet  with  a  Solid  Core. 
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Figure  3c.  Evaporation  Data  for  a  Mixture  Particle  Initially  Composed  of  25  %  Syloid 
72  and  75  %  Octyl  Alcohol  and  Calculated  Evaporation  Rate  of  an  Octyl  Alcohol  Droplet 
with  a  Solid  Core. 
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Figure  4a.  Comparison  of  Evaporation  Data  for  Mixture  Particle  Depicted  in  Fig.3a 
and  Shrinking  Core  Evaporation  Model  Calculations  for  Three  Values  of  the  Effective 
Diffusion  Coefficient  DE.  "Best  fit"  DE=0.015  cmVsec. 
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Figure  4c.  Comparison  of  Evaporation  Data  for  Mixture  Particle  Depicted  in  Fig. 3c 
and  Shrinking  Core  Evaporation  Model  Calculations  for  Three  Values  of  the  Effective 
Diffusion  Coefficient  DE.  "Best  fit"  DE=0.025  cmVsec. 
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Notations 


c°  =  octyl  alcohol  saturation  vapor  concentration  at  temperature  T  (moles/cm^) 
Cex=  octyl  alcohol  vapor  concentration  external  to  particle  (moles/cm^) 

Cj  =  octyl  alcohol  vapor  concentration  inside  porous  shell  (moles/cm®) 
d  =  particle  diameter  (cm) 

Dg  =  octyl  alcohol  gas  phase  diffusion  coefficient  (cmVsec) 

De  =  porous  shell  diffusion  coefficient  (cmVsec) 

Dp  =  pore  diffusion  coefficient  (cmVsec) 

m  =  particle  mass  (gm) 

m^  =  initial  particle  mass  (gm) 

m,  =  silica  mass  (gm) 

mi  =  octyl  alcohol  liquid  mass  (gm) 

Mw  =  octyl  alcohol  molecular  weight  (gm/mole) 

N  =  a  constant  (moles/sec) 

Pd  =  octyl  alcohol  partial  pressure  at  droplet  surface  (ergs/cm^) 

Po»  =  octyl  alcohol  partial  pressure  at  infinity  (ergs/cm*) 
r  =  radial  coordinate  (cm) 
r^  =  pore  radius  (cm) 

R  =  radius  of  solid  particle  (cm) 

Rj  =  gas  constant  (ergs/mole- °K) 
t  =  time  (sec) 

T  =  temperature  (°K) 

V,  =  octyl  alcohol  volume  (cm*) 
v,=  octyl  alcohol  molar  volume  (cm*) 

V,  =  silica  volume  (cm*) 

V  =  levitation  voltage  (volts) 

V,,  =  initial  levitation  voltage  (volts) 

Vd  =  dry  particle  levitation  voltage  (volts) 

V,  =  levitation  voltage  of  octyl  alcohol-saturated  particle  (volts) 

Greek  Letters 

6  =  free  volume  of  particle  (cm*/gm) 

K  =  solid  tortuosity 

17  =  solid  surface  porosity 

p,=  silica  density  (gm/cm*) 

Pi=  octyl  alcohol  density  (gm/cm*) 

T  =  shell  thickness  (cm) 
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